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We theoretically investigate intermodulation distortion in high-Tc superconductors. We 
study the effect of nonmagnetic impurities on the real and imaginary parts of nonlinear 
conductivity. The nonlinear conductivity is proportional to the inverse of temperature owing 
to the dependence of the damping effect on energy, which arises from the phase shift deviating 
from the unitary limit. It is shown that the final-states interaction makes the real part 
predominant over the imaginary part. These effects have not been included in previous 
theories based on the two-fluid model, enabling a consistent explanation for the experiments 
with the rf and dc fields. 
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The high-temperature superconductors are attractive for use in microwave circuits because 
of their low surface resistance as compared to those of normal metals.^ This low-loss property 
is disturbed by the nonlinearity of its response to external fields. This nonlinearity means that 
the system is unstable with respect to changes in the input power. This effect makes super- 
conductors unsuitable for practical applications. On the other hand, the nonlinear response 
is useful to investigate the intrinsic properties of superconductivity. It has been predicted 
that the nonlinear Meissner effect (NLME) shows a peculiar behavior in unconventional su- 
perconductors.^ This prediction is summarized in the following two points. One is that the 
nonlinear correction to the magnetic field penetration depth (A) is proportional to the inverse 
of temperature (T). Then, the divergence at low temperatures yields a nonanalytic response. 
The other is that the nonlinear correction takes different values depending on the direction of 
the external field. These can be evidence of the existence of nodes in superconductors. 

Several experiments have been conducted on this effect; these experiments yield different 
results depending on the measurement methods. An experiment that measures the dependence 
of A on the magnetic field yielded a result that is inconsistent with the theoretical prediction.^ 
(Neither the low-temperature upturn nor the angle dependence is observed.) Intermodulation 
distortion (IMD) is theoretically supposed to reflect the presence of the NLME.^ An experi- 
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ment on the IMD seemingly shows a result consistent with the theoretical prediction.^ (Only 
the low-temperature upturn is observed. The angle dependence has not been investigated. In 
this sense, this experiment is inadequate to be considered as evidence for the NLME.) 

In this paper, we microscopically develop the theory of nonlinear microwave response 
and consider how this contradiction arises. In the IMD experiment, the power is measured, 



which is expressed as Pimd |Ai?s + iAXs\'^ Here, AR^ = Aai / {2a2) to j-i / cr2 and AXg = 
— A(T2/(2cj2)\/'^/f^/o"2 for CJ2 >> O"!. (Rg, Xs, and a = ai — \a2 are the surface resistance, 
surface reactance, and conductivity, respectively, and A implies the nonlinear correction.) 
Previous theories on IMD assume the validity of the two-fluid model in addition to that of 
the theory of Yip and Sauls; in these theories, only Acj2 is considered.'^ (We show that this 
assumption does not necessarily hold.) This is suflicient when the response to the nonlinear 
dc field'^ is considered. In the case of the IMD, however, there is a contribution from Aai in 
general. 

In the linear response, it is known that a2» (Ti holds. On the other hand, the relationship 
between Aai and Acj2 is not known. Therefore, we calculate both real and imaginary parts of 
the nonlinear conductivity to determine which quantity is predominant. We have to specify 
a dissipation mechanism in order to estimate Adi, though this is not the case for Ac72- The 
NLME conies into question at the low-temperature region, where the 1/T-upturn is supposed 
to occur. Therefore, we mainly consider the effect of nonmagnetic impurities on the nonlinear 
microwave response. This is because with regard to dissipation, the impurity scattering effect 
is dominant at low temperatures and the electron-electron correlation is dominant near T^-^ In 
this sense, we do not consider the type of correlation effect that functions as the enhancement 
factor and can be effective in the response to a static external fleld.*^ The absence of the 
NLME under a nonlinear dc field can be explained by taking this effect into account. An 
explanation to the above contradictory behavior can be provided by combining this effect 
with the invalidity of the two-fluid model discussed here. 

We consider isotropic impurity scattering. The self-energy with the self-consistent t-matrix 
approximation is 

^'^'>- cot'6-G^{er 
Here, Ti = ni/7rN{0) (n^ and N{0) are the impurity density and the density of states at the 
Fermi level in the normal state, respectively) and Gq^c) = Tr G^^/(27rA^(0)) with the 
Green function 

rR - 1 (' + ^' (2) 

(e = e— EQ^(e).) The nonlinear response function (third order) can be expressed as follows. (The 
vertex correction is given by the functional derivative of the self-energy by the one-particle 
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Green function as that in the conserving approximation/'' which is also derived from Keldysh's 
method on the nonequihbrium state. "'^^) K^^^ui, 102,^3) = / dei^e \u)i,U!j,U!k ). (The 

conductivity is expressed as Aa oc K^^^/co.) X^^j^ ;-] means the sum of ah permutations 
{i,j, k} = {1, 2, 3} and lo = coi + L02 + i^3- 

(f f \;^RAAA I f ^AAAA 

— yjel — Je2)9ele2e3e4 + Jel9eU2e3e4 

_i_T^r.r f ShRRR IRRR , iRRR r)RR tRRR \ ff f \tRAR t^RR fRRR 

"r-'-r[— /e4t'ie4ele2-'^e2,e4''e2e3e4 "r "e3e4cl^cl,c3'''clc2c3 J ~ \ Je3 — Je4j'ie3£4el^cl,c3'''clc2c3 

af \i,ARA n^A iAAA , f jiAAA ri^A iAAA , iAAA nAA iAAA \] 
el — /e2j'ie4ei£2-'-^e2,e4"'e2e3e4 "r /ell"e4ele2-'-^e2,e4"e2e3e4 "r '^e3e4el-^el,e3'''ele2e3/J 

-L-rvr tf f \hARR n^A tRRA ff f \tARR n^A tRAA 
+ -'-r[-(,/e3 - /e4jfte4ele2-^e2,e4"'e2e3e4 " Ue2 - /e3 j'ie4ele2-^e2,e4"e2e3e4 

ff f \tAAR T^RA IRRA (f f \tAAR j^RA fRAA i 

- Ue2 - /e3j'ie3e4el-'^el,e3"eU2e3 ~ l/el " /e2j'ie3e4el-'-'el,e3"ele2€3j- 

(3) 

Here, ^,7,™ = -;^Gfi:.-.Gf2'.-fcGf33,.,Gf,^,, Z.™ ^ G7.-.Gf2^.-..Gf3^., is 
the quasiparticle velocity, = tanh(e/2T), ei = e, €2 = e — uji, = e — uji — uj2, £4 = e — a; , 
and LJ = aji+a;2+i^3- and i^^^ are vertex corrections; they are given afterward. The first 
and second traces represent the variations of the density of states and self-energy under the 
external field, respectively. The third trace implies the vertex correction, which represents the 
final-states interaction. The reasons for the invalidity of the application of the two-fluid model 
to the nonlinear response are as follows. (1) It is based on the assumption that the damping 
effect is independent of energy. (2) It includes only the nonlinear response of the density 
of states (the dependence of the damping effect on the external field and the final-states 
interaction are omitted). Therefore, we investigate these two aspects. The diagrams of the 
nonlinear response are shown in Fig. 1. Fig. 1(a) shows the diagram of the linear response; 
the nonlinear corrections are shown in (b) and (c). Fig. 1(b) and (c) show the nonlinear 
response with the variation of the density of states and vertex correction, respectively. In the 
linear response, vertex correction does not exist in the case of isotropic impurity scattering. 

First, we consider the nonlinear response arising from the variation of the density of states; 
its response function can be expressed as follows. 



CO 



^/d#>(^)i. (4) 



de 3 \ de^ J 7, 

■^2^wR 



Here, 7^ = -lmSo(e), n^^ = /pg V^e/^e^ - A|, and Vb,u,«; = '^,'"l,vl. (We substitute = 
Aocos20 and take Aq as the unit of energy in the following numerical calculations.) If 7^ is 
independent of energy, we have the same result as that when the two-fluid model is used. 

(3) (3) 

The temperature dependences of ReKj-,Qg/uj and lmKj-,Qg are shown in Fig. 2. As the 
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(b) 




(d) 



Fig. 1. (b) and (c) Representative diagrams for K^^\ The solid and wavy lines denote the one-particle 
Green function and the external field, respectively. The shaded rectangles denote vertex correction, 
(a) Linear response diagram, (d) Diagram of the impurity scattering effect with the self-consistent 
t-matrix approximation. 
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Fig. 2. Temperature dependences of Rei^j^Q^/wjoj^o and Imi^jj^^ for various values of the phase 
shift 5 and Ti = O.OOL 



,(3) 



phase shift S deviates from unitary scattering {S = 7t/2), ReKj^Qg/col^^^o attains larger values 
and becomes proportional to the inverse of temperature. On the other hand, lmK)-^Qg does not 
show a clear l/T-divergence, but it is cut off at low temperatures. (The graph oi S = 0.4757r 
is seemingly divergent, but this is also verified to be cut off by comparing with that of smaller 
S or 1/T.) This behavior of K^^q^ can be explained by the dependence of the damping rate 
on energy. In previous theories on the nonlinear response in the Meissner state, the l/T- 
divergence is supposed to arise from the derivative of the density of states, which is cut off at 
low temperatures by the impurity scattering.^^ (In clean systems, Ked'^ne/de^ oc S{e) because 
Rerig oc |e|.) If the damping rate takes a constant value, the result shown in Fig. 2 cannot be 
explained. The energy dependence of the damping rate 7^ is shown in Fig. 3. As S deviates 
from 7r/2, 7^ decreases at around e c± 0. Then, Red^nf^ / de^ increases, but it is cut off at low 
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Fig. 3. Energy dependence of 7^ with the same parameters as those used in Fig. 2. 



energies because 70 7^ 0. Therefore, the different dependences on temperature arise in Fig. 2. 
'ReK^^Q^/oj shows a 1/T-divergence owing to the energy-dependent damping effect. Imif^^^ 
is cut off at low temperatures reflecting the energy dependence of Red'^nf^/de^. 

The nonlinear correction to K^^^ resulting from the variation of the self-energy is shown 
as a diagram similar to that in Fig. 1(c). The four-point vertex is expressed as D^^ = 
niT^^/{l — UiT^'^'mdn^^ / de). This term is small as compared to the vertex correction -D^^, 
which is verified by a numerical calculation. Therefore, we omit this term. 

Next, we consider the contribution of vertex correction to K^^\ which is written as 

K% = J de^y ^ M^i + iu:N2)D^^^^{N^ + iu;feiV2). (6) 

Here, Ni = Re{dnl^ / de) h„ N2 = [Ni/-fe + lm{d'^n^^^/de^)]/{2-f,), and D^^^^^ = {lo - coi + 
2ij^)/[K{uj — u>i)Yien^^/j^]. (The term with D^^ does not exist in the case of a nonlinear 
dc field.) The way in which the vertex correction depends on frequency originates from the 
identity 



'7rAr(0) 



(7) 



(here, f,^ = {-cotSh - G§JttN{0))-^ and 



which is similar to the 



1 ' 

identity discussed in the localization problem. ^"^ 

In the numerical calculation of the two-tone IMD we substitute cji = 0^2 = + 
and c<j3 = —u — 2Auj and then maintain 101^2.3/^ as constant for uj ^ 0. The contributions 



from the vertex correction, YicKy^/oj\^^Q and ImKyl^, are shown in Fig. 4. At 6 = 7t/2, both 



(3) 



ReKyQ/ioluj^o and ImKy^, decrease as the temperature decreases. As the phase shift deviates 
from tt/2, they show an upturn as 1/T increases. These behaviors are explained by the energy 
dependence of the damping rate and its effect on the density of states. Both KeKy^/io\i_j^o 
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(3) 

and ImKyQ are independent of the phase shifts at high temperatures. This means that the 
impurity scattering effect is less dependent on the phase shifts in this temperature region, 
as shown in the high-energy part of Fig. 3. The dependences of K^^^ on phase shifts appear 
in the low-temperature region. The expression of K'^^^ indicates that ReK^]~,/Lo\^^o and 
ImKy^ are proportional to and 7°, respectively. This distinguishes the behaviors of the 
real and imaginary parts of K^^\ ReKy(~,/uj\ui-*o shows an almost 1/T-divergence, but ImK^l 
is roughly proportional to 1/ Vt. The absence of the cut off at low temperatures in liaKyQ 
(unlike the case of ImK^Qg) originates from the energy dependence of the density of states 
(rig^) in the vertex correction D^^. 

In Figs. 2 and 4, we can see that the real part of K^^^ shows a l/T-divergence at some values 
of phase shifts. On the other hand, ImK^^^ does not show such a behavior. We should clarify 
which of Aai and Aa2 is predominant in order to specify the origin of the low-temperature 
upturn in the IMD power. To see this, we evaluate the following ratio. jolKeK^^^ /lo) /ImK^^\ 
which is equivalent to (7o/a;)Acri/A(T2, is shown in Fig. 5. In the hydrodynamic regime, which 
is the premise of our calculation, 70 is greater than u. Therefore, Aa2 is always predominant 
over A(Ti if (7o/a;)Acri/A(T2 < 1 holds. On the other hand, there is a possibility of Acri > Aa2 
in the case of (7o/a;)Acri/A(T2 > 1, depending on the value of 70/a;. As shown in Fig. 5, if 

(3) 

we consider only Kj-^Qg, Aa2 > Aai holds in the same way as the two-fluid model. When we 

(3) (3) 

take account of Ky(j, Aai can predominate Acr2- As can be seen from Figs. 2 and 4, ImKy^y 
takes values of the same order as ImK^Qg. On the other hand, ReKyj^/u takes values that 
are 100 times greater than ReK^j^g/uj. This difference originates from the following fact. It 
can be shown that the term D^"^ oc l/Au arises in the real part of K^jy (this term is cut 
off by the nonlocal effect mentioned below), but this term is canceled out in the imaginary 
part. Therefore, Acri can possibly predominate A(T2; as a result, the l/T-divergence can be 
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Fig. 5. Temperature dependences of ^oCReK^^ /i^)/lTxiK^g^ (here, S = DOS or DOS + VC) with 
various values of phase shifts and the same parameters as those used in Fig. 4. 

originated from Aui. This yields a solution for the contradiction between the experiments 
with the nonlinear rf and dc fields, which is not resolved when the two-fluid model is used. 

Here, we mention some issues that are not discussed above. Strictly speaking, in the case 
of 5 7^ 7r/2, the self-energy should be written as the matrix S^(e) = So^(e)fo + S3^(e)f3. 
(S3^(e) = — rjCot5/[cot^(5 — Go^(e)^] and fo is a unit matrix.) In this paper, we present the 
formula with (e) — because it yields an intricate expression of K^^\ and this gives 
almost the same numerical results as those for (e) ^ 0. We present the numerical results 
of K^^^ calculated by using the original expressions ij^^ie) ^ 0). The diamagnetic terms that 
include the factor dvk/dkrs are omitted. This is because out of the two branches, only the 
gap-full branch remains in the vertex correction D^^; these branches arise from the matrix 
structure in the superconducting state. The nonlocal effect is not considered here. This effect 
also broadens the singular behavior of the derivative of the density of states in the same way as 
that by the impurity scattering effect.^' However, the thickness of the film used in the IMD 
experiments^ is nearly 4000 A, which is almost 100 times thinner than that of the experiment 
with the nonlinear dc field.^ This is almost the same order of magnitude as A. Therefore, we 
omitted this effect here. (The numerical calculation of the current distribution with various 
values of A is given in ref. 15.) We show only the numerical results in which the impurity 
concentration Fj was fixed. This is because our argument on K^^^ can be similarly discussed 
when Fj is varied. The different points are that the phase shift at which the l/T-divergence 
appears depends on Fj and the absolute value of K^^^ varies with Fj. 

In our theory, whether Pjmd oc l/T^ or not depends on the value of the phase shift, which 
is not known so far. As for the phase shift deviating from the unitary limit, however, there 
are several discussions related to the low-temperature thermal conductivity that suggests 
neither unitary nor Born limits. ^^'^^ With regard to the comparison between the real and 
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imaginary parts of the nonlinear conductivity, there is an experimental suggestion that ARg is 
predominant over AXs,^^ though 1/T-diver gence is not expected to exist in their temperature 
range. One of the possible experiments that can verify our theory is the third harmonic 
generation. When oji = uj2 = 1^3, the contribution from vertex correction to Rei^(^) is reduced 
to the same order as Rei^jj^^. Therefore, it is expected that C72 > (Ti holds and the l/T- 
divergence is cut off at low temperatures. 

In this paper, we derived the general formalism of the nonlinear microwave conductivity 
under the influence of nonmagnetic impurities. We evaluated this formula by varying the value 
of the impurity scattering phase shift. As the phase shift deviates from the unitary limit, the 
nonlinear response shows a 1/T-divergence owing to the dependence of the damping rate on 
energy. This is one of differences from previous theories where the 1/T-divergence originates 
from the second derivative of the density of states. The predominance of the resistive part over 
the reactive part arises when the vertex correction is included. This term is not included in the 
two-fluid model. Therefore, the upturn of the IMD power at low temperatures can originate 
in the resistive part. This upturn does not need to be accompanied with 1/T-divergence in 
the reactive part; this is a possible explanation to the seemingly contradictory results between 
the static and microwave experiments. 

The numerical computation in this study was carried out at the Yukawa Institute Com- 
puter Facility. 
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